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ABSTRACT 



Data representing the image of an original picture having 
high resolution are first acquired. The data are then pro- 
cessed in accordance with a blurring function, which can be 
a randomly distributed point spread function. Processing can 
be carried out either in object space or in spatial-frequency 
space. In object space a multi-dimensional convolution is 
required In the spatial frequency domain, processing is 
performed by a pointwise product of data representing the 
Fourier transform of the original image and the Fourier 
transform of the blurring function. The resulting blurred 
image is thereby coded such that the original scene is no 
longer discernible, but can be recovered provided the blur- 
ring function is known. Additionally, this blurred image can 
be shown to contain the detail of the original but can be 
stored or transmitted as data which are compressed by 
compression ratios which can range from 2 to 100. The 
compressed data can be restored to provide an image with 
detail like the original by inverse filtering using data repre- 
senting the Fourier transform of the blurring function. 

31 Claims, 4 Drawing Sheets 
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IMAGE DATA CODING AND COMPRESSION 
SYSTEM UTILIZING CONTROLLED 
BLURRING 



DESCRIPTION 

The U.S. Government has rights in this invention pursuant 
to contract DAAL-03-92-G-0147 between the Department 
of the Army and the University of Rochester. 10 

This application is a continuation-in-part of our applica- 
tion Ser. No. 08/095,328 filed Jul. 21, 1993. 

The present invention relates to image data coding and 
compression systems wherein the data representing an origi- 
nal image or picture can be transmitted or stored as a coded 
and compressed representation which renders the original 
scene unrecognizable and also requires far fewer pixels 
and/or bits than the original uncoded image. 

The invention is especially suitable for use in compress- 1Q 
ing pictorial images into a representation having a compres- 
sion ratio which can be selected to range from about 2 to 
100. Compression may be achieved through either of two 
means: (1) by subsampling the data thereby reducing the 
number of pixels in the compressed representation, or (2) ^ 
reducing the quantity of data used to represent the image by 
reducing the average number of bits-per-pixel in the repre- 
senting data. The compression ratio is defined as the ratio of 
the number of pixels in the original to the number of pixels 
in the representing compressed data for the case of sub- 30 
sampled data or the ratio of the original quantity of data to 
the quantity of data needed in the compressed representation 
for the case of a bit reduction system. The compressed data 
represent a coding and compression of the original image in 
accordance with a selected random code, which unless the 35 
code is known, ii not recoverable; this .coded form being 
particularly advantageous for storage, transmission, and 
privacy. 

This invention is especially useful for coding a picture 
prior to transmission or storage. For example the picture is 40 
greatly blurred, physically it takes on the appearance of a 
photograph recorded through turbulence. The controlled 
blurring greatly reduces the visual perception of edges and 
high frequency content is no longer discernible thereby 
rendering the original image unrecognizable. One may wish 45 
to code the picture for privacy reasons alone and this 
application is not dependent on also using the image com- 
pression aspects of this invention. 

The present invention is based on the discovery that 
controlled blurring of an image provides information which 50 
contains all of the detail of the original and yet allows fast, 
accurate recovery of the original. One can show that the 
blurring causes a radical restructuring of the data in the 
spatial transform domain, greatly increasing amplitudes at 
low frequencies relative to high frequencies. In the case of 55 
a subsampling system, the number of pixels required to store 
or transmit the blurred image is then much smaller than for 
the original image. For the bit reduction system, the reduc- 
tion in high frequency energy increases the compression 
ratio due to the reduction in bits required to code the high 60 
frequency information. In other words, the blurred image 
can be stored or transmitted as a compressed representation 
of the original image, without sacrificing the detail of the 
original image even though the blurred image appears 
degraded to an extent that the original is not recognizable 65 
from the blurred version. The amount of blurring and 
compression depends upon the blurring function. The blur- 



2 

ring function is statistically a distribution characterized by a 
probability density function. Hence an infinite array of 
specific realizations can be generated for coding or blurring. 
One can generate a specific realization using a pseudo- 
random process generator. A particularly desirable function 
is a randomly distributed set of impulses yielding in trans- 
form space a speckle pattern or in general a random-walk 
transfer function. The function can have a uniform distri- 
bution or an exponential distribution or it can be generated 
using a fractal formula. The distribution may be Gaussian. It 
can be shown that the compression ratios can theoretically 
exceed 100. Compression ratios from 2 to 100 have been 
obtained in typical cases of interest. 

Another desirable function for blurring is an impulse 
response consisting of the sum of two terms: an impulse at 
the origin and a Gaussian spot which is broad enough to span 
several pixels. For even more secure coding a final illustra- 
tive example is a point-spread-function having the following 
terms: a delta-function or impulse at the origin and a set of 
randomly arrayed Gaussian spots. 

Compression and coding of image data has involved 
techniques operative on the data itself such as transforma- 
tion coding using various types of Fourier transforms such 
as the discrete-cosine-transform (DCT). See M. Rabbani and 
S. Daly, Proc. Soc. Photo-Opt. Instru. Eng., 1071, 246-256 
(1989) and U.S. Pat. No. 4,703,349 issued OcL 27, 1987 to 
J. G. Bernstein. Other techniques have involved filtering of 
signals obtained from the image thereby reducing the fre- 
quency content and thereafter restoring the frequency con- 
tent to retrieve the image. See, U.S. Pat. No. 4,777,620 
issued Oct. 11, 1988 to Schimoni, et al. The literature 
relating to image reconstruction does not discuss blurring 
using the multiple impulse point spread function nor the 
random- walk transfer function and subsequent partial image 
reconstruction. See chapter entitled 'Image Enhancement 
and Restoration - by-B. R. Frieden in the treatise; Picture 
Processing and Digital Filtering, ed. by T S. Huang, 
Springer-Verlag, Berlin 1979, and particularly Table 5.1 
thereof. It has also been proposed to send a pilot signal 
subject to the same encoding transformations as the image 
representing signal through the encoding process. The pilot 
signal is separately detected and used to compensate for 
losses in picture information which are incidental to the 
encoding process. See U.S. Pat. No. 4,807,029 issued Feb. 
21, 1989 to Tenaka. None of the processes discussed in the 
above-referenced patents nor any other processes for encod- 
ing have involved blurring of the image or using the blurred 
image to provide compressed data. It has not heretofore been 
appreciated that the blurred image contains all the detail of 
the original and that the blurring function provides the key 
or code by means of which the original can be restored from 
the blurred version without significant loss of detail or 
resolution. 

It may be appreciated that an image formed by superim- 
posed exposures on a photosensitive receptor of the original 
at different displacements results in a blurred image on the 
photosensitive receptor. Observation of this image makes it 
apparent that the detail, represented by high frequency 
content of the original, is not visually perceived in the 
blurred picture. In frequency space the high frequency 
information is attenuated and the area containing non- 
negligible information represented by the blurred picture is 
smaller than the corresponding area for the original image. 
Accordingly, the original image can be represented by 
compressed data. The distribution of the displacements is 
preferably a random distribution. The locations in which the 
multiple exposures are made as well as the number of 
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exposures controls the blurring. Knowledge of the displace- 
ments (the blurring code) may be used in an inverse process 
to recover the original image. The process of blurring and * 
recovery is preferably carried out with the aid of a digital 
computer or computers. Compressed data representing the 5 
blurred image may be transmitted or stored thereby reducing 
the time for transmission or memory space for storage over 
what would be required for transmission or storage of the 
original. 

Accordingly, it is the principal object of this invention to 10 
provide an improved system for coding and compressing an 
image and particularly data representing that image whereby 
high compression ratios can be obtained and the original 
image accurately recovered from the compressed represen- 
tation thereof. 15 

It is a further object of the present invention to provide an 
improved system whereby data may be encoded into a 
compressed form so as to provide compression ratios which 
can be selected in accordance with the encoding function. 

T M^ a stjll fiirthpr object of the prymt inYHitirm t rt 20 
provide an improved system for d ata comnrfj^jn" " f imagft 
representing data using blurred versions of the ima ge 
thereby maintaining fidfllity pi the image while requiri ng 
fewer pixels than in the original image , and enabling resto- 
ration or decoding only if the code for blurring is known. 25 

It is a still further object of the present invention to 
provide an improved system for data encoding and possible 
compression of image representing data which utilizes con- 
trolled blurring in accordance with a blurring code which 
renders the original scene unrecognizable in the encoded 
image, can be decoded only with knowledge of the blurring 
function^ and can readily be changed from time to time so as 
to suit different security requirements. 

Briefly described, a system for representing an original 
image by^compressed data from which a restored version of 
the original can be produced with approximately the same 
detail as the original comprises means for spatially coding 
the original to provide a blurred image. Means are then 
provided for spatially decoding the blurred image in accor- 
dance with the blurring code to provide the restored version 
of the original image with like resolution and detail, 

The foregoing and other objects of the invention, as well 
as personally preferred embodiments thereof, will become 
more apparent from a reading of the following description in 45 
connection with the accompanying drawings in which: 

FIG. 1 is a block diagram illustrating a system for data 
compression and restoration utilizing a blurring function to 
process a pictorial image with a controlled blur, compressing 
the blurred picture by means of either subsampling or 50 
scanning the blurred image or utilizing less bits to represent 
the blurred data, transmitting or storing of the compressed 
data, and decompressing and deblurring to restore the origi- 
nal image; 

FIG. 2 is a block diagram or flow chart illustrating the data 
handling processes performed by the components of the 
system shown in FIG. 1; 

FIG. 3 is a block diagram illustrating one form of the 
compressor shown in FIG. 1 utilizing a subsampling system; 

FIG. 4 is a block diagram illustrating one form of the 
compressor shown in FIG. 1 utilizing a bit reduction system; 
and 

FIG. 5 is a flow chart illustrating the processes used in 
computing the blurred image which are carried out in the 55 
digital data processor of FIG. 1. 

Referring to FIG. 1, an original pictorial image denoted 
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35 



40 



55 



60 



functionally by p(x,y) in x and y orthogonal coordinates, is 
acquired by image acquisition using a digitizer 10. The 
digitizer may have a video camera which may be a CCD 
array which provides a bit map of nxn pixels representing 
the image. This map is stored in memory in a digital 
processor 12 (a computer). The blurring function b(x,y) is 
inputted to the computer from a random pattern generator 
14. The computer 12 itself may be programmed to generate 
the pattern. The pattern is two-dimensional and is, in a 
preferred form of the invention, a set of N impulses at 
positions (x n ,y„) chosen from a random process having a 
probability density function f(x„) f(y„). These points may be 
from a uniform density function, Gaussian function (in 
general, an exponential function), or any of a wide variety of 
probability functions known to those versed in statistics. 
These functions represent random codes. IVpical codes for 
different compression ratios are shown in Table I which is 
presented and discussed below. In effect, these random 
patterns produce, in transform space, speckle patterns which 
constitute a random-walk type of transfer function. The 
number of points, N, as well as the distribution of the N 
points determines the compression ratio. The distribution 
may be selected depending upon the desired compression 
ratio using 25-250 points. A compression ratio of pixels in 
the original image to pixels in the blurred image (the 
compressed data) can range from 2 to 100. These compres- 
sion ratios can be obtained with 25 points spread over a 
64x64 pixel array. For the case of a subsampling system, the 
compression ratio obtained then is determined by the sam- 
pling density at which the blurred picture is sampled in a 
sampler type data compressor 16a, in FIG. 3. 

The processor 12 is programmed to convolve the original 
image "data p(x,y) with the blurring function, b(x,y) to 
produce data representing the blurred image,, g fc (x,y). It can 
be shown, and is shown hereinafter, that the blurred picture 
requires- less volume -in frequency space no represent 
adequately the image than the unblurred picture and can be 
represented by compressed data. The blurred picture is in the 
form of a bit map which may be compressed using any 
standard bit reducing algorithm or it may be sampled by the 
sampler so as to provide compressed data. For the latter, the 
number of pixels is equal to nxn divided by the compression 
ratio, n being the number of pixels along each dimension of 
the original. 

The compressor 16 in the case of a subsampling system 
may be a coarse resolution facsimile machine of the type 
which is commercially available. Then the image which is 
scanned is a pictorial image which may be formed on a 
screen or xerographically as hard copy. A bit map in storage 
in the computer 12 may also be sampled coarsely to provide 
the compressed data. Alternatively for the case of a bit 
reduction system, the compressor, as shown at 16a, FIG. 4, 
may use other compression methods operative on digital 
data representing the blurred image, g fr (x,y). Such tech- 
niques include the DCT discussed in the Rabbani, et al. 
article referenced above and predictor methods such as 
differential pulse code modulation (DPCM). In either case 
the compressed data is applied to a transmitter 18 for 
transmission over a link, such as a radio or telephone 
channel. Alternatively the data may be stored and forwarded 
on the link. The compressed data may also be stored either 
optically or electronically on various recording media (e.g., 
magnetic tape or discs), and later recovered. 

In the case of transmitted, compressed data, a receiver 20 
acquires the transmitted data and stores it as the blurred 
picture data g b {x t y). A digital data processor 22 decom- 
presses and deblurs the data utilizing a blurring function 



08/25/2004, EAST Version: 1.4.1 



5,453,844 



20 



25 



30 



generator 24, which may be part of the processor 22. Since 
the blurring function b(x,y) is known, deconvolution or 
inverse filtering obtains the restored image, p(x, y) which can 
be displayed on a graphics display 26. 

In the case of a subsampling system, the compressor 16 
reduces the number of pixels to be transmitted or stored. 
This sampler as shown at 16a in FIG. 3, utilizes a display of 
the blurred image g fr (x,y), which may be a video display or 
a photographic or xerographic copy of the blurred image. 
The display is indicated generally at 28 in FIG. 3. The 
display is scanned by a coarse scanner 30, such as a low cost 
fax machine, which provides the compressed data to a 
modem 32 (modulator which may have a demodulating 
function) which transmits the data over the link. 

The processes performed in the computer 12 and in the 
processor 22 are shown in FIGS. 2-4. The original sharp 
picture p(x,y) is convolved with the blurring function b(x,y) 
by a two-dimensional (2D) convolution to provide the 
blurred picture data g fc (x,y). This blurred picture is stored in 
two dimensions in an nxn pixel array. In the subsampling 
system the sampler then transmits and/or stores compressed 
data by coarse scanning the number of pixels transmitted, 
n*xn*. The number of pixels in the product n*xn* may be, 
for example, Vi to Viooth the number of pixels in the nxn 
array. It is this compressed data that is transmitted over the 
link and/or stored. This compressed data is received or read 
from storage. Then a 2D Fourier transform is performed on 
the compressed data resulting in Fourier transform space, 
G b (f x , The Fourier transform of the blurring function 
B(f x , f y ) is obtained by another Fourier transform process. 
These^Fourier transform processes may both be fast Fourier 
transforms (FFT). 

The blurred picture transform G^f* i y ) is divided by the 35 
blurring transform BCf^, fp and the quotient is applied to a 
Fourier inverter also operative in two dimensions. The result 
is the original sharp picture p(x,y) which may be displayed 
or printed. It can be shown that the zeros of B(f x f y ) are 
isolated and thus there is a neglible likelihood that a zero can 40 
occur to disrupt the inversion process. 

The convolution process utilizes a 2D-FFT to obtain the 
Fourier transform of the image, ?(f x ,f y ). Another 2D Fourier 
transform is used to provide the blurring function transform 
B^fy). The random process generator produces N pairs 
x„,y„ which can be mathematically represented by equ ation 
(5) which is illustrated in FIG. 5 at the input to the FFT 
process which produces an output Btf^fy). 

A multiplier produces the product of P(f x ,fp and Btf^y. 
This product is subject to a two-dimensional inverse FFT to 
produce the blurred image g & (x,y) which is then sampled to 
provide compressed data for storage or transmission. 

The following presentation shows that the blurred version 
of the original image requires significantly fewer pixels for 
transmission or storage then the original hence providing 
data compression ratios for example from 2 to 100. The 
restoration of the original image from the compressed data 
uses the blurring function in an inverse filtering process as 
discussed above. 

Consider an image p(x,y) that has a frame size X by Y and 
a pixel size Ay=Ax. Based on the Whittaker-Shannon sam- 
pling theorem, (see, e.g., "Introduction to Fourier Optics", 
Joseph W. Goodman (McGraw-Hill Book Co., San Fran- 
cisco, 1968) Section 2.3 p. 21 and problem 2.12 page 28.), 
the number of pixels required to characterize this picture is 
given by Np as 



45 



XY 
AxAy 



(1) 



In the frequency space using Nyquist sampling concepts, 
the image p(x,y) occupies in spatial-frequency (or Fourier 
transform) space denoted by A that is given by the following: 



A=F Jt F y 

in which 



l 

Ar 



(2) 



(3) 



1 

by 



The size of the Fourier transform space, in spatial fre- 
quency coordinates, is F x by F y . Now, substitution of Eq. (3) 
into Eq. (1) yields the number of pixels expressed as the 
product of the area of the picture in (x,y) space and the area 
in spatial-frequency space, viz. 



(4) 



From the foregoing discussion the space-bandwidth prod- 
uct of a compressed version of the original image after 
blurring may be derived. Consider the blurred image as 
defined by the point-spread function (psf) that consists of N 
randomly distributed points (x n ,y„). The original image is 
characterized in terms of its intensity by p(x,y) which is a 
real-valued function of real (x,y) and it is also non-negative. 

The blur function b(x,y) is given by in which N is an 
integer typically from 2 to 250. For purposes of 



AT 
«=1 



x n )6(y-y 0 ) 



(5) 



this presentation N=25. The numbers (x rt ,y„) are obtained 
from using the random-process-generator. Statistically x n 
and y n are chosen from distributions characterized by a 
probability density function, f(x„) f(y n ). The numbers x„ and 
y n are taken as uncorrected and both as having the same 
probability density function. 

The blurred picture is denoted by g^x.y) and is given by 
the expression 



50 



■jj: 



(6) 



P^~x\y-y , )b{x\yWH' 



In Eq. (6) this integral form is often denoted by an asterisk, 
as in Eq. (7), which shows Eq. (6) as a two-dimensional 
convolution: 
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g fc (*,y)=p(x,y)*b(x,y) 



(7) 



The blurred image g fr (x,y) contains sufficient information 
to enable recovery of the original image, p(x,y,),; very 
accurately. Moreover, the number of pixels required to 
characterize g^(x,y) is considerably less than that required to 
characterize p(x,y); thus providing significant data compres- 
sion for storage or transmission of data representing the 
original image and from which it can be recovered accu- 
rately. This compression will be apparent from the reduction 
in the space-bandwidth product of g^x.y) relative to p(x,y). 

Consider the two-dimensional Fourier transforms of both 
sides of Eq. (7); the result is 
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c»(f*y=P(f^B(f,,y 



(8) 



where capital letters denote the Fourier transforms of the 
corresponding functions. The blurring function B(fjpf y ) acts 
as a low-pass filter in the spatial-frequency domain. Con- 
sider the function b(x,y); it has a Fourier transform B(f„fy) 
defined by the following integral: 



(9) 



10 



Substitution of Eq. (5) into Eq. (9) yields the following form 
for the transform of the blurring function: 



(10) 



Using statistical analysis, we calculate the expected value of 



8 



<BB*> is small and so there is little need for storing this 
data. Thus the band-limiting envelope (in frequency space) 
controls the compression. The number of points N deter- 
mines the amount of high frequency reduction since high 
frequencies are reduced by a factor of 1/N relative to low 
frequencies. 

Table 1 contains the bandlimiting envelope that controls 
the compression, <BB*>, for two important practical 
choices of the random code: (1) Uniform density and (2) 
Gaussian density. Two different criteria are useful to 
describe the controlled blurring. These arc tabulated in 
columns: (a) using bandwidth to the first null or to 1/e and 
(b) more conservative limit defined as 1/N 2 of the second 
term of Eq. (14). 



TABLE 1 




Typical Blurring Estimates 






Density Function 


Compression Envelope 


Compression Criteria 


for Random Code 


<BB*> 


(a) 


(b) 


1) Uniform 

2) Gaussian 


N + N(N + l)sinc a (sfjsinc 2 (sf y ), 
s = mAx . 
N + N(N+ l)e- (lBoiS(fl2 ^ 
a = qAx 


On) 2 
(2jcq) 2 


(rcm) 2 /|N(N + 1)) 
(2rcq) 2 /ln[N(N+ 1)| 



B(f,,f y ), denoted by the ensemble average bracket notation: 

B<f^)>=NF(yP(f y ) (Ll) 

In Eq. (11), the function F(f x ) is the characteristic function 
(the Fourier transform) corresponding to the probability 
density f(x J. This notation is as follows: 
The value of the spectral density is computed using Eq. (8). 
This appropriate form (where (f^fy) are omitted for 



= J dxrj{x n ) 



<e Hlnr »*»> = F{Q 
simplification) is 

<G t G/>=PP*<BB*> 



(13) 



<BB*>=N+N(N-lW<f x )\ z \FVhd y)l 2 



<BB*>=N 2 when f x =f y =Q 
<BB*>=N. as f„f,->~ 



35 



A generalized controlled blurring function may also be 
used. This function is described by 



40 



45 



Again, the band-limiting property is seen by exarnination of 
the term <BB*>. From Eq. (10) <BB*> can be expressed as: 



(M) 



where F(f y ) is the characteristic function of the probability 
density f (y„), derived similarly toF(Q as shown in Eq. (12). 
Critical values of <BB*> using Eq. (14) are as follows: 



55 



(15) 
(16) 



In other words, in spatial frequency space, Eqs. (13) to (16) 
show the compression of the original data PP*. The shape of 
the envelope <BB*> controls the distribution of energy in 
the blurred image and thus controls the amount of compres- 
sion achievable. At large f x and f y , the envelope function 



N 



71=1 



Z s(x-x R ,y-y n ) 



(17) 



(18) 



50 



where s(x,y) is a sampling function such as the gaussian or 
a triangle function. This function provides greater smooth- 
ing and hence better compression and better reduction of 
perceived detail in the coded image. For illustration, s(x,y) 
can be written as follows: 



3{x,y)=c 



(19) 



in which w is a constant which is arbitrarily chosen, but in 
a magnitude range suitable for computer operations. 

Another blurring function can be generated by using an 
impulse response of the following form: 



60 



b 2 {x,y)=Cd(x)b{y)+b l {x.y). 



(20) 



In frequency space the equivalent multiplicative transfcr 
function is given by taking the two dimensional Fourier 
transform of b 2 (x,y); the result for B 2 (f x f y ) is given by 



65 BtfJ^C+BtfJJ 

Where only one exponential is employed i.e., 



(21) 
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_ *y (22) 

fcC*.y) = C5tx)8(y) + e 

The shape of the bjfoy) function in frequency space is a flat 5 
pedestal capped by the exponential function for low spatial 
frequencies. The blurred image is very smooth and contains 
only small high frequency components. These high fre- 
quency components being about an order of magnitude 
lower than the low spatial frequency components. Trie 
compression of such blurred images is particularly efficient 
using Discrete Cosine Transform (DCT) or differential pulse 
code modulation (DPCM) methods. 

From the foregoing description, it will be apparent that 
there has been provided an improved data compression 
system utilizing controlled blurring with a blurring function 15 
which is a random code thereby also providing security of 
transmission and/or storage. It will be appreciated that 
different blurring codes, which can be changed from time to 
time may be used depending upon the data compression 
ratio and degree of secrecy required. Variations and modi- 20 
fications in the herein-described system, within the scope of 
the invention, will undoubtedly suggest themselves to those 
skilled in the art. Accordingly, the foregoing description 
should be taken as illustrative and not in a limiting sense. 

We claim: 25 

1. A system for compressing and decompressing image 
data signals representing an original image characterized by 
p(x,y), where x and y are the spatial coordinates of pixels of 
the image, which system comprises: 

(a) means to convolve said image data signals with coded 30 
random pattern data signals generated according to a 
blurring function 



(b) said recovering means includes means for providing, 
in response to the compressed data signals representing 
G*(f*r»f y ) and data signals representing B(f,,L), data ' 
representing the quotient thereof, and means for pro- 
viding the inverse Fourier transform of said quotient to 
provide said data signals representing said original 
image in restored form. 

5. The system according to claim 1 further comprising 
means operative according to a random process for obtain- 
ing x n and y„. 

6. The system according to claim 5 wherein said random 
process operative means is a generator providing a generally 
random point spread function of N points in an area in x and 
y containing said image. 

7. The system according to claim 6 wherein N is in a range 
from 2 to 250. 

8. The system according to claim 1 wherein said means to 
compress said convolved data signals further comprises 
either means for subsampling said signals after convolving 
with the blurring function or means for reducing the number 
of bits representing pixels at each of said spatial coordinates 
x &y. 

9. The system according to claim 8 wherein said data 
compressing means operative on pixels of said blurred 
image is selected from the group consisting of means for 
subsampling said pixels and means for reducing the number 
of bits of data signals representing said blurred pixels. 

10. The system according to claim 1 wherein the blurring 
function is another blurring function which is described by 



b\{x,y) = s(x,y)* I 8(x-*n)S(y-yn) 
n=l 

N 



b(x,y) = I :${x-x n )h<y-y n ) 



35 



where N is an integer greater than one and x rt and y n are 
the positions of delta-function impulses, x„ and y n 



being chosen from distributions characterized by prob- 
ability density functions, to provide convolved data 
signals representing a blurred image; 

(b) means to compress said convolved data signals to 
permit efficient storage or transmission of the informa- 45 
tion in said original image; and 

(c) means to recover said original image from said con- 
volved compressed data signals with a function which 
is the inverse of said compressing means and said 
coded blurring means to provide data signals represent- 50 
ing said original image in restored form. 

2. The system according to claim 1 wherein said signal 
convolving means, said compressing means, and said recov- 
ering means are operative in object space, said convolving 
means providing the convolution, g fr (x,y>=p(x,y)*b(x,y). 55 

3. The system according to claim 1 wherein said signal 
convolving means, said compressing means, and said recov- 
ering means are operative in spatial frequency space. 

4. The system according to claim 3 wherein 

(a) said signal convolving means and compressing means 60 
comprises means for providing the Fourier transform of 
the data signals representing p(x,y) and b(x,y) which in 
said spatial frequency space are pCf^fy) and B^fy), 
respectively; means for multiplying P(f K) f y ) and B(f x ,f y ) 
to provide the product thereof, G^.fy); and means for 65 
providing said compressed convolved data signals in 
response to and representing G b (f xt f y ); and 



where s(x,y) is a sampling function including Gaussian and 
triangle functions. 
11. The system according to claim 10 wherein s(x,y) is 



40 



jfey)=e 



in which w is a constant. 

12. The system according to claim 10 wherein one of the 
impulses is taken at the origin i.e. 5(x)8(y) and includes the 
point (Xj^Ojj^ti). 

13. The system according to claim 10 wherein an impulse 
function at the origin C8(x)8(y) is added to the generalized 
controlled blurring function bj(x,y) to give b 2 (x,y) defined 
by: 

b&y^WxXlyKbfry). 

14. The system according to claim 13 wherein the blurring 
function is: 



W*,y) = CS(*)5(y) + c * 

15. The system according to claim 1 wherein an impulse 
at the origin C8(x)8(y) is added to define the blurring 
function. 

16. The system according to claim 1 further comprising 
means to store or to transmit said compressed convolved 
data signals. 

17. A system for representing an original image having a 
resolution characterized by a first array of picture elements 
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(pixels) with a second array of pixels less in number than the 
pixels in said first array and from which said second array a 
restored version of said original image can be reproduced 
with approximately the same detail as said original image, 
said system comprising means for spatially coding said 5 
original image according to a spatial code to provide a 
blurred image, and means for spatially decoding said blurred 
image responsive to said spatial code to provide said 
restored version of said original image. 

18. The system according to claim 17 further comprising 10 
means for translating said blurred image into said second 
array, and said spatial decoding means including means 
responsive to said spatial code which codes said image for 
translating said second array into an array of pixels generally 
corresponding to said first array. 15 

19. The system according to claim 18 wherein said spatial 
code is a generally random point spread function of N points 
where N is an integer greater than one. 

20. The system according to claim 19 wherein N is in a 
range from 2 to 250. 20 

21. The system according to claim 18 wherein said 
translating means comprises means for scanning said blurred 
image to provide a pattern of pixels forming said second 
array. 

22. The system according to claim 18 wherein said 25 
translating means comprises means for sampling said 
blurred image to provide data signals representing said 
second array. 

23. The system according to claim 17 wherein said spatial 
coding means includes means operative upon the Fourier 30 
transform of said original image in accordance with a 
selected speckle pattern to provide said blurred image, said 
spatial decoding means includes means operative on said 
blurred image in accordance with said selected speckle 
pattern. 35 

24. The "system according to claim 23 wherein said 
speckle pattern is represented by the Fourier transform of a 
plurality of randomly distributed points, N in number where 
N is an integer greater than one. 

25. The system according to claim 9 further comprising 40 
means to store or to transmit said second array of pixels. 

26. A method for representing an original image having a 
resolution characterized by a first array of picture elements 
(pixels) with a second array of pixels less in number than the 
pixels in said first array and from, which said second array a 45 
restored version of said original image can be reproduced 
with approximately the same detail as said original image 
comprising the steps of: 



spatially coding said original image according to a spatial 

code to provide a blurred image, and 
spatially decoding said blurred image in response to said 

spatial code to provide said restored version of said 

original image. 

27. A method according to claim 26 further comprising 
translating said blurred image into said second array, and 
where said spatially decoding step includes translating said 
second array into an array of pixels generally corresponding 
to said first array. 

28. A method according to claim 26 wherein said spatial 
code is a generally random point spread function of N points 
where N is an integer greater than one. 

29. A method according to claim 26 wherein said step of 
spatially coding according to said spatial code includes 
operating upon the Fourier transform of said original image 
in accordance with a selected speckle pattern to provide said 
blurred image, and said step of spatially decoding in 
response to said spatial code includes operating on said 
blurred image in accordance with said selected speckle 
pattern. 

30. A method according to claim 26 whereby said spatial 
code is a blurring function represented by 



«=1 



yn) 



where x and y are spatial coordinates of pixels, N is an 
integer greater than one and x„ and y„ are the positions of 
delta-function impulses, x n and y n being chosen from dis- 
tributions characterised by probablity density functions, to 
provide convolved data signals representing a blurred 
image. 

31. A method according to claim 30 wherein said blurring 
function is described by 



b)(x,y) - jfcy)* I 8(x - x„)6(y - y„) 
«=1 

N 

= Z >s&-x n ,y-yn) 
n=L 

where s(x,y) is a sampling function including Gaussian and 
triangle functions. 

* * * * * 



50 



55 



60 



65 



08/25/2004, EAST Version: 1.4.1 



